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Abstract 

An algorithm for the reduction of one-loop n-point tensor integrals to basic integrals is 
proposed. We transform tensor integrals to scalar integrals with shifted dimension [jl], 
and reduce these by recurrence relations to integrals in generic dimension §. Also the 
integration-by-parts method [Q] is used to reduce indices (powers of scalar propagators) 
of the scalar diagrams. The obtained recurrence relations for one-loop integrals are 
explicitly evaluated for 5- and 6-point functions. In the latter case the corresponding 
Gram determinant vanishes identically for d = 4, which greatly simplifies the application 
of the recurrence relations. 



1 On leave of absence from JINR, 141980 Dubna (Moscow Region), Russian Federation. 
2 Supported by BMBF under contract PH/05-7BI92P 9 



1 Introduction 



Radiative corrections contain not only most essential information about quantum properties of 
a quantum field theory but moreover, their knowledge is indispensable for the interpretation of 
precision experiments. Higher order calculations in general have been performed for processes 
with at most four external legs. With increasing energy the observation of multi particle events 
is becoming more and more important and at least one-loop diagrams must be calculated 
for these cases as well. One of the important examples is the W^-pair production reaction 
e + e~ — > W + W~ — > 4 fermions (7), being experimentally investigated at LEP2, which allows 
us to accurately determine mass, width and couplings of the IV-boson. At future high- 
energy and high-luminosity e + e~ linear colliders radiative corrections will become even more 
important for a detailed understanding of Standard Model processes. Thus their computation 
will be crucial for a precise investigation of particle properties and the expected discovery of 
new physics. At present, a full 0(a) calculation for a four-fermion production process is not 
available, mainly because, such a calculation is complicated for various reasons [^, |J. One 
of the problems is addressed in the present article and concerns the calculation of one-loop 
diagrams with special emphasis on processes with five and six external legs. 

The evaluation of one-loop Feynman graph amplitudes has a long history - [19| and at the 



present time many different methods and approaches do exist. However, the most frequently 
used Passarino-Veltman JTIJ approach is rather difficult to use in calculating diagrams with 



five, six or more external legs. The tensor structure of such diagrams is rather complicated. 
To obtain the coefficients in the tensor decomposition of multi-leg integrals one needs to solve 
algebraic equations of high order which often are not tractable even with the help of modern 
computer algebra systems. Additional complications occur in the case when some kinematic 
Gram determinants are zero. 

In the present paper we propose an approach which simplifies the evaluation of one-loop 
diagrams in an essential way. It allows one to evaluate multi-leg integrals efficiently without 
solving systems of algebraic equations. Integrals with five and six external legs are worked out 
explicitly. Those with seven and more external legs need special investigation because there 
are different types of Gram determinants vanishing and their consideration is postponed for 
this reason. The above cases with more than four external legs allow particular simplifications: 
for n=5 due to the property (|5l|) and the property of the recurrence relation (B7I) and for n > 5 
due to properties of Gram determinants. Our scheme allows also to evaluate diagrams with 
n < 4, but here we cannot use the above mentioned properties. A simplification is achieved 
in this case since solving of linear equations is avoided and the algorithm is implemented in a 



FORM |20| program for these cases. 



2 Recurrence relations for n-point one-loop integrals 



First we consider scalar one-loop integrals depending on n — 1 independent external momenta: 

^ = J^fljj. (i) 

J 71 j=i c j 







where 



Cj — (q — Pj) 2 — m 2 + ie for j < n and c n = q 2 — + ie. 



(2) 



The corresponding diagram and the convention for the momenta are given in Fig.l. 

9-W ' w,- F3 
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Fig. 1: One-loop diagram with n external legs. 
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can be written as a combination of scalar integrals with shifted space-time dimension multiplied 
by tensor structures in terms of external momenta and the metric tensor. This was shown in 
Ref.|l| for the one- loop case and in Ref. for an arbitrary case. We shall use the relation 
proposed in || 



(4) 



where T is a tensor operator, dj = and d + is the operator shifting the value of the 

space-time dimension of the integral by two units: d + I^ = I^ d+2 \ On the right-hand side 
of (^) it is assumed that the invariants q have arbitrary masses and only after differentiation 
with respect to m 2 these are set to their concrete values. 

To derive an explicit expression for the tensor operator T^ 1 „^ r {{p s ) ) {dj} ) d + ) we introduce 
an auxiliary vector a and write the tensor structure of the integrand as 



?w • • ■ % 
Next we transform the integral 
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into the a-parametric representation by means of 
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(5) 



(6) 



(7) 



and perform the <i-dimensional Gaussian integration 



d d kexp i(Ak 2 + 2(pk)) 
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The final result is: 
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where 



(8) 
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Q(a,{p s },a) = + 5I(ap fc )a fc - -a 2 



fc=i 



(10) 
(11) 



and Q(a, {p s }) is the usual a- independent Q-form of the graph. From the representation (||) 
together with (ED it is straightforward to work out that 



^,.^(fe},{^} ) d + ) = in-^exp 
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This representation is particularly well suited and effective for a computer implementation of 
the tensor integrals (0). 



2.1 Integrals with non zero Gram determinants 

The purpose of this Sect, is to develop an algorithm for reducing the above mentioned scalar 
integrals to standard integrals in generic dimension d = 4 — 2s. 

Recurrence relations which reduce the index of the j-th line without changing the space-time 
dimension are obtained by the integration-by-parts method H: 



2A r ^.j+/W = £(l + ^ 
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( dA n 
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(13) 



where is the Kronecker delta symbol, the operators ^ etc. shift the indices z/,- 



and 
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Taking derivatives of A n one should consider all Y^- as independent variables and for j > k 
assume d/dYj k = d/dY k j. After taking derivatives one should set 



Yij = ~(pi - Pj) 2 + m 2 + m 2 



(14) 



where Pi,Pj are external momenta flowing through i— ,j-th lines, respectively, and rrij is the 
mass of the j-th line (p n = 0). At this stage the external momenta are not yet restricted to 
dimension 4. We will specify later when this property is used. 



A recurrence relation for reducing dimension and index of the j-th line is obtained from 



(<9,A n ) + £(<9 j( 9 fe A„)k- 

k=l 



± n i 



where dj = djdm 2 - and 
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(16) 



We may also reduce the space-time dimension of by means of: 



i=l 



2A n + J2(dkA n )k- 



k=l 
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(17) 



Equations (]IB|), ( ]T5] ) and ( [T7D are our starting point. Some simplifications of these equations 
can still be obtained. In particular we give in the following a compact representation for the 
mass-derivatives of A n . First of all we mention the following useful relation between G n _i and 



A n : 



3=1 



-Gn-l- (18) 

The basic object for the purpose of expressing the derivatives in (0), QI5]) and (^) turns 



out to be the "modified Cayley determinant" of the diagram with internal lines 1 . . . n |10] , 
namely 
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labeling elements . . . n. "Signed minors" 



J1J2 • • • 
Mk 2 • • ., 



(20) 



will be labeled by the rows ji,j2, ■ ■ ■ and columns kx, k 2 , ■ ■ ■ excluded from () n [110], |2"2"fl . E.g. 
we have 

A„ = (°) (21) 



Further relations are 



dA n /dY jk = (2-5 



djA n 



did k A r 
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(22) 
(23) 
(24) 



and for p n = we have 

Gn-i = 2 () n - (25) 

One obvious advantage of the above relations is that they can be easily used for numerical 
evaluation. Recursion ( |13) then reads 
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we write it in the most convenient form for further evaluation as 
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1 + 



k-i)|^-e A!/ 



(28) 

Here the 'deviations' [y k — 1) of the indices from 1 are explicitly separated so that all indices 
z/fc = 1 do not contribute in the second sum in curly brackets on the r.h.s. Finally the double 
sum can be completely reduced to a single sum by means of |], || 



3=1 



(29) 



This relation reduces simultaneously indices and dimension, which is what one wants in gen- 
eral. It is not possible, however, to introduce it directly into fl28|) since in (|28| ) we explicitly 
have to separate the term i = k. Therefore, further details of how to apply recurrence relation 
(13) can be given only in the case of explicit examples (see Sect. 3). In our notation recurrence 
relation fll5|) now reads 



and recurrence relation ( |T7|) 
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(31) 



Relations (|2g) (including @), (0) and (0) now replace QT3|),([T^) and (0). For z/, = 1 (|5D|) 



and (|3lD correspond to eqs. (20) and (18) of the first paper of jl5| . In Sect. 3 we demonstrate 



how to reduce tensor integrals to scalar ones, to which these recurrence relations are then 
applied. 

Furthermore the following general observation is needed in what follows. Applying the recur- 
sion relations, contractions of the n-th line occur. In this case we encounter integrals with 
p n 7^ 0. In order to apply our recursion relations for this case, we must properly define G„_i. 
We can either shift all momenta by p n or we can use from the very beginning the definition of 
G n _i in terms of flT9| ) with p n ^ 0. In this manner the recurrence relations remain unchanged. 

If our recurrence relations are to be implemented in terms of a computer-algebra program, 
the following steps are recommended. At first, by using (Q) we reduce tensor integrals to 
scalar ones with shifted dimension (see also below). In a second step we apply ( |l5l , |3~0"l) and 



then (|T7],[n]). Both relations produce the same factor G n ^\ in the denominator and therefore 



one expects some similarity of the obtained expressions. Finally, relations (|T3|j28| ) are applied, 
bringing all integrals to a set of master integrals in the generic dimension with powers 1 of the 
scalar propagators. 

The properties mentioned in the introduction leading to simplifications for n > 5, which will be 
explained in more details below, are difficult to be implemented in computer algebra systems. 
Therefore we prefer specific representations, e.g. to avoid Gram determinants Q In particular, 
to obtain results as compact as possible, we will do a careful analysis of the recurrence relations 
and different properties of Gram determinants in Sect. 3. 



2.2 Integrals with zero kinematic determinants 

When one or both of the determinants G n _i and A n are equal to zero the reduction procedure 
must be modified. In such cases it is possible to express integrals with n lines as a combination 
of integrals with n — 1 lines. 

First let us consider G n _i = 0. This is the case for n-point functions with n > 6 for external 
4-dimensional vectors since then the order of () ra is n + 1 but its rank is 6. In this case the 
recurrence relations flHUP and ([JT]) cannot be directly used. Prior to the reduction of d one can 
remove one of the lines of the diagram by using a relation which follows from (|15|) : 



J « ) = _^(^n) k - J W, i (32) 
k=i l a i ZA «J 

Here we keep the original form in terms of derivatives w.r.t. the masses since in some computer 
algebra systems these derivatives may be easier to calculate than determinants. 

By repeated application of this relation one of the lines will be contracted. The procedure 
must be repeated until one obtains integrals with non-vanishing determinant G. After that 
the reduction of the space-time dimension can be done for integrals with the smaller number 



3 Strictly speaking, here we only mean G n _ 
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of lines. For G n _i = yet another relation can be obtained from flI7D: 



^ = -E^k-e (33) 
k=i ZLAn 

In fact both decompositions are equivalent since for G n _i = 0, cancelling the common factor 
in A n and <9fcA n , the following relation holds 

<9 fc A n = d k djA n 

2A n ~ d, t An ' 1 j 

We will see in the next Sect, that indeed further important consequences follow from G n _i = 0. 



Zeros of A n may occur in special domains of the phase space for n=6. In this case the reduction 
procedure will be as follows. First, applying the relation 

n n 

(d- 5> - 1)G„_!/W = J2(d k A n )k-lt 2 \ (35) 

i=l k=l 

which follows from (pj]) , the integral with n propagators must be reduced to a sum of integrals 
with n — 1 propagators but in lower space-time dimension. As in the case with G n _i = the 
procedure of contracting lines can be repeated until we obtain integrals with non- vanishing A n . 
After that integrals with space-time dimension d < 4 can be expressed in terms of integrals in 



the generic dimension d by using (|29|) and indices can be reduced to 1 by using relations ([13]) 
and (B3). 



Explicit calculation ^ shows that A n = for n-point functions with n > 7, again for external 
4-dimensional vectors - and in particular for arbitrary masses. Then, due to ( ^3|) also (^j = 
for n > 7. These properties eliminate quite many terms in relations fl28|), (|30|) and (|31|) 0. This 



is the reason why the investigation of n > 7 is postponed at present. Nevertheless, we write 
the recurrence relations in what follows for arbitrary n. In particular for external vectors not 
of dimension d = 4 the Gram determinants do not vanish as described and for these cases 
anyway our relations hold in the given form for arbitrary n. 



3 Explicit recursions 

In this section we give details of how to calculate in particular 5- and 6- point functions since 
these are of greatest actuality for present day experiments. Our goal is to present the most 
compact formulae for direct applications. Throughout we use the notation 

d d q . . 

(36) 



7T 



d/2 



4 This has been performed for n = 7 and 8 in terms of a FORM program, using component representations 
of fourvectors. 

5 See also the discussion after (pQ). 



s 



and for the scalar integrals we introduce the explicit notation 



j[d+Y,stu... __ I Y[ C -( 1 + S ri+S rj +5 rk + ...-8rs-SrtSru-.-) (37) 

■> r=l 

where [d+] 1 = 4 + 21 — 2e . Observe that ( number of entries s,t,u, . . . ^ i,j,k, . . .)+p — n and 
equal upper and lower indices 'cancel'. The index p specifies the 'actual' number of external 
legs. We consider now Feynman diagrams in the 'generic' dimension d = 4 — 2s. For scalar 
n-point integrals in the generic dimension we use the notation I n . 



For the reduction of the tensor integrals to scalar ones with shifted dimension, we wrote a 
FORM [20f program which applies @ to integrals of rank 1, 2 and 3. After inspection the 
result reads for arbitrary n and powers of the scalar propagators equal 1 (the latter is the 
most frequent case in the electroweak Standard Model in the Feynman gauge; otherwise, as 
mentioned above, we can reduce higher indices by recurrence relations): 



^ n ^ 

r=l 

n-1 

= Y.V^nf (38) 
1=1 

r d n 

C = / sVIK 1 

J r =l 

= ZW - V/j^ (39) 
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i,j,k=l 
n-1 



i n— i 

+2 E orrf + + , (40) 
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where = 1 + 6ij and = 1 + 6ij + 6ik + 6jk — SijSik^jk is the number of equal indices 
among i,j, k which can be written in this symmetric manner. 



We first consider now n — 5, using recursion relation (|30|) we have: 

,. T [d+] 3 _ 5 r[i+l 2 ,f 6)r[«i+P,» _ 1 I X _|_ X (AU 

"ijk h.ijk 7\ + 2^ T\~ 2 P,ij ' ^ - J- + ^ifc + (41 J 

Us 5=1 Us 

where in the second term on the r.h.s. we have introduced the index p. In general p — 4, 
but for s = i,j we have p = 5 (recall that equal upper and lower indices cancel). Let us now 



o 



consider the first term on the r.h.s. of ( f4"T| ) 
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(42) 
(43) 
(44) 



where we indicated two more steps of the reduction. Like in flilD we have again introduced 
the index p in the second term on the r.h.s. of ([42]): for s = i, p = 5. At this point it is also 
worth pointing out that 

f^ij ^iji ^ijJfc' ^ij^ijk' (45) 

Now consider the second term on the r.h.s. of fl4"T|): 
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(46) 



Here the index g in the second term on the r.h.s. is q = 3 in general, except when t = i, in 
which case g = 4 (see (0))- F° r 9 = 3 we have: 
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l 3,i 
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(47) 



For s — i,j in flip ) and also in the second part of (fE]) we have integrals of the type 



T [d+r 

1 5,i 



Wj T 
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5 E () 



0/5 r[dH 
J- a 



(48) 
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and finally, applying recursion (p~7|) , we have 



r[d+] 



r [d+],s 



8=1 



5 ts 



(0s\ 5 /0s\ 



(:) 



t=i 



0, 



rst 
'3 



d-4 
1 



d- 3 



(49) 



(50) 



We observe that fl4*9|) is unpleasant in the sense that the expression in the square bracket 
is zero for d = 4 and the overall factor is — ^, i.e. we need to expand the 5- and 4- point 
functions up to order e in order to get the finite part. In fact, however, we will show that 
/I cancels. To demonstrate this for the 5-point function, we need first of all to express the 
g^ v tensor in terms of the external vectors. Assuming that the external vectors Pi, ■ ■ ■ ,Pa are 
4-dimensional and independent (i.e. no collinearities occur), we can write ( see e.g. 



flU 



2£ 

ij=i 







(51) 



1 n 



Inserting this into (BDI), we get 



E m ■ n^ffS - E Wry ■ ( 52 ) 

By inspection of (H^) we see that the second term on the r.h.s. for s — i exactly contains the 
to cancel the in (|52|). In a similar way with the help of ( [51"D we can rewrite the 
second sum in fllPP , to get 



Eu v A 
PiPjPk 



i,j,k=l 



S ki 0. J 0." . 



(53) 



To pick out the /^"^ we use ([48]). Inserting ( |41~D into (40), the J^" 1 "' contributions come from 



the second term in (f42|) (for s = i) and the second term of (|4l|) (for s = i, j). Further we need 
the property ([45]). 



Finally a remark is in order concerning the fmiteness of the 5-point function: the only infinities 
occurring in the above decomposition are coming from the I2 's in (0). The | terms in these 
2-point functions are independent of masses and momenta, however. Thus according to (see 

0) 



Ef ) = > i = l,-'-n (54) 



for n = 3 these infinite terms cancel. With these considerations the tensor 5-point functions 
can be completely reduced to scalar 4-, 3- and 2-point functions in generic dimension. We 
want to point out that use has been made only of relations ( |30"D and (pi]) . We will see be- 
low that (f43D should be replaced by another relation containing no Gram determinant, which 
might simplify the numerics in some kinematic domains. At the end the tensor integrals un- 
der consideration are finite (after cancellation of the above mentioned terms of order | ) and 
therefore we can put e = 0, which in particular applies to the scalar 5-point function 0, [1~0|, P~5[| . 

For n = 6 and d = 4 the situation is completely different due to the fact that ()$ = as was 
discussed in Sect. 2.3 . Therefore (|15|) and ( |I7|) both reduce to (^), i.e. instead of those two 



recursion relations we have only one. Explicitly it reads, introducing the signed minors of the 
modified Cayley determinant 

, 6 ( R ) 

= E TST^S!: 1 ' = any value 0,..., 6. (55) 

We see that d55|) does not allow to reduce the dimension and it can only be used to reduce 
indices! If p = 5 on the r.h.s. of (|55|), which is the case for r 7^ i.., then the reduction can be 
continued as for the 5-point function above. If, however, p = 6 on the r.h.s. of (|5q ) in case of 
r = i.., then again ( p5|) has to be applied until all Iq L are eliminated. 



So far only relations fl3!]| ) and ( |3"T| ) for the Gram determinants G n have been applied. Now 
let us investigate (p8|) . It does not contain Gram determinants G n . A priori it does not allow 

1 1 



to reduce the dimension. In the following, however, we will explicitly show that ( |2"9"| ) can be 
used to perform one sub-summation in (p8j) and that it is this relation which finally provides 
a very efficient possibility to reduce also the dimension in the case of the 6-point function. 
Also some further useful information about the reduction of 5-point functions is obtained in 
this manner. In the following we investigate explicitly the integrals occurring in (|38l),([39|) and 
Hit). 



For arbitrary n (pE|) yields for the integral in (|38| ) 
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and using (|29| ) once we obtain : 



(56) 



AX 



[n + l-(d+2i)]H e ]i + E 



r [d+](i-i) )7 . 
J n-1 



Similarly, for the integral in (|39|) we obtain from (28) for arbitrary 
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(59) 



It is interesting to note that in (E7|) for n = 5 and I = 1 as well as in ([5^) for n = 6 and 1 = 2, 
the numerical square brackets evaluate to 4 — d = 2e. This means a great simplification for 
e = and in particular (0) ought to replace (|3|) for n = 5 ( a little algebra shows that indeed 
for n = 5 (^) and (i3) are identical including the £-part). Moreover this is exactly what is 



needed in (^) for n = 6 and 1 = 2, the integral needed in Q3"9] ) for n = 6. 
For the integral with three indices occurring in ([40"D we finally have 
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Equations (|5TD,(|59D and (|60D are valid for any n and they represent in a way an optimal form 
of the recursion for the considered integrals. For n > 7 and d = 4, as mentioned in Sect. 2.3 



i o 



already, these relations are still valid but reduce considerably. Since this case is at present of 
minor physical interest, it will be discussed separately. We only mention that, e.g., using (|60|) 

for i = j = k and n = 7 we obtain a relation for l\ d ^ in terms of ijf^ ' r , summed over r, etc. 

The integrals with highest dimension in the above relations have coefficients (^j , (^j and 

(fy , respectively. In (|38|), (|39|) and fl4"0|) these are multiplied with pf, pj and p£ and summed 
over i,j,k. Due to 

eW°) =o. n > Q ( 61 ) 

all these contributions vanish. To prove (|6T|), we project it on all pj, i = 1 ... n — l(p n = 0) in 
d=4 dimensions, assuming that at least four p/s are linearly independent. Thus we have to 
show 

n fo\ 

Eftft ■ =0 , t = l...n-l. (62) 

3=1 V/ n 

First of all we write the scalar products in the form (] 



so that 



reads 



PiPj 2 Yin Yjn ^nn} 



o 



(63) 



(64) 



^ ] {Yij Y n j {Yin Y n n)} 
j=l V/ n 

Now the term (Yi n — Y nn ) is independent of j and the summation over j can therefore be 
performed: 

n /Q N 



E 

i=i 



).- 


()» = 


o, 


n > 6. 


1) 








Ya 


**2 


^3 


y. 

• • • 1 in 


Y n 


Vl2 


^13 


Yin 


Y X2 


^22 


Y23 


Y-2 n 


Y tl 


*i2 


Y i3 


Y 

■ ■ ■ 1 in 


Yin 


^2n 


Ys n 


Y 

• • • 1 nn 



(65) 



,0, 



(66) 



Subtracting the 0-th line of this determinant from the i-th one, the latter takes the form 
(1, 0, ... 0). Expanding finally the determinant w.r.t. the first column yields the r.h.s. of fl66"D. 
Thus this determinant is independent of i, which finally proves fl64j) , and hence (|6l"r) . 

Further we also see that ^ in (|59"D cancels against the g^ u term in (|39|) (I = 2) due to 

2 "=1 /0 z' 



n > 6. 



(67) 



6 This essential decomposition, to our knowledge, was first used in (l(J and is the basis of the proof in jl5| 
as well. 
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As above we show this again by proving that all projections of this tensor on any pair of 
vectors Pk,Pi{k, I — 1 . . -n — 1) vanish, i.e. 

Expressing all scalar products again in terms of the Yik, this is equivalent to 

2J Kfe — — {Xkn — Y nn )] \Yj£ ~ Yj n — (Xtn ~ ^nn)] I n . ] 
i,j=l \ U -Vn 

= [Y H - Y kn - {Y ln - Y nn )} Q (69) 
The contribution on the l.h.s. of ( |55| ) with K's independent of z, j is proportional to 

sfiP.-s©-().- ft 

The 'linear' terms, according to ( f27|) and (|66|), read 

i*£M--ffl pi) 



i=l 



and are independent of k, i.e. the contributions of this type cancel. Finally, the 'quadratic' 
terms are of the form 

HHS3-£ y '-M!).-M!5.' (72) 

which proves ( p9|) and thus (|67|). 

We point out that (|67|) is a representation of the 5^ tensor in close analogy to (0), only that 
in this case the n — 1 vectors (n > 6) are linearly dependent in <i = 4 dimensions. 

Finally we investigate the simplification of fl40| ) due to (0), where we concentrate on the 
6-point function, n = 6. Let us first consider the I n d p integrals. These are given by (157|) for 



/ = 2. First of all we again observe that their first part can be dropped due to (61). Their 



remaining parts for n = 6 are only of the I5 J type, which we would like to cancel. This 



concerns in particular the complete second part in (f40|) . It is easy to see that due to (|67|) the 
g^ x and g uX terms cancel against the two terms in the square bracket of (^) by summing in 
Of) over i, k and j, k , respectively. To show the cancellation of the remaining /| d s is a bit 



tricky. At first we observe that the last term in (|60|) contains the following term (most 
easily seen from fl42|) ): 

T [d+] 2 ,r _ Qr) T [d+],r 
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Thus the sum of all remaining terms is given by 



Then we observe that for n = 6 the tensor can be written in the following form 



!^ = 2 EW. r = 1...6. (75) 

For r < 5 this equation immediately follows from fl5"T|), just by scratching one of the 5 mo- 
menta and assuming that always four of them are linearly independent. For r = 6, (|75|) can 
be proven in the same way as fl67]). Thus, summing over i,j = 1 ... 5 in flT4"|), the cancellation 
of all /J^ -1- ' integrals is also shown for the tensor integrals of rank 3, a result which has been 



obtained previously in fll5| . What remains in ( p0[) is finally again only the last sum with the 
exclusion of the contribution. The ,r integrals in (|60| ) can be calculated, e.g., by 
means of (|42|) , without the term s = i in the sum on its r.h.s. or from (|59| ) again dropping 
on the r.h.s. . 



As a final remark we point out that the only Gram determinants occurring are those in the 
tensor integrals of rank 3, coming from the reduction in (|42|), where instead of (|43"D now (|57|) 
is to be used. Thus only inverse Gram determinants of the first power appear. Moreover the 
remaining integrals are all reduced to standard 4- and 5-point functions for n = 6. 

Conclusion 

We have presented an algorithm for the calculation of n-point Feynman diagrams, applicable 
for any tensorial structure and gave a detailed specification for 5- and 6-point functions. 
Dependent on the kinematic situation different recurrence relations may be used. A typical 
example is the equivalence of (f43"D and (|57|) for n = 5. All recursion relations are valid for 
arbitrary dimension. In the main part of the paper we have concentrated on the experimentally 
most relevant cases, i.e. the tensor 5- and 6- point functions, and considered d = 4 for the 
external momenta in order to simplify the results. In this case a regulator mass is needed for 
infrared divergences. We restricted ourselves to tensors of rank 3 []. For the 5-point function, 
due to the above equivalence, there occur inverse Gram determinants of only second order; 
for the 6-point function, surprisingly, only of first order, n-point functions with n > 7 are not 
considered explicitly since due to the drastic reduction of the recurrence relations for d = 4 
these cases need a separate investigation. 

Acknowledgment 

We are grateful to A. Davydychev and S. Dittmaier for useful discussions and to A. D. also 
for carefully reading the manuscript. O.V. Tarasov gratefully acknowledges financial support 
from the BMBF. 

7 For more than three integration momenta in the numerator the extension of our procedure is 
straightforward . 
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Appendix 



In this appendix we present for the 6-point function the 'effective' contributions to the integrals 
( p8|) - (0), i.e. leaving out all those contributions which we have shown to cancel. First of 
all no contributions come from the metric tensor in ([59]) and (|4Tj|). The integral in Q5^| ) is 
'effectively' given by 



1 6,i 



1 



r=l 



.Or, 



1 5 



1 y- (or) 

r=l 



E 



'0r N 



jrs 
M 5 



(76) 



where we have used ([)7|) and expressed the 5-point function in generic dimension in terms 
of 4-point functions. No Gram determinant is involved in this case. 



The remaining integral in ([39]) is 'effectively 



n--I [d+]2 



a? 



r [<*+], r 



= i 



oV =1 



6 ^, 



Or 
Or J 6 



E 



v 0sr . 



TTS 

1 A ■ 



(77) 



Here (|59p is used and the integral 1^ is expressed in terms of (0), where the square bracket 
of the first contribution evaluates to 2e and is dropped. The result is very similar to (|76|), 
i.e. it is a double sum over 4-point functions and no Gram determinant involved. The latter 
property is due to fl57|) , which is used instead of (|43|). 



The integral in 



is a bit more complicated: 



■n 1 T^ s 



1 



6 

E 

r=l 



Or 
Or 
Qr 



w..j[ d +]V 



'Or 



r[d+] 2 ,r 



L 5.i 



+ E 

s=l 



'Ojr 
v 0sr , 



r [d+],rs 
'4,i 



(78) 



Here ( pf^) is used together with fl59|) for the reduction of ,r . In the square bracket of the 
first part of ( |59|) again e = is taken. Reducing if^ ,r we see no possibility to avoid the 
inverse Gram determinant and therefore propose to use fl48|) , i.e. 'effectively' (dropping 



L 5,i 



1 6 

E 



with (e = 0) 



r [d+],rs 



\rsj ( 



[Qrs, 



sr 



r [d+],rs 



jrs \ * 
J 4 



jrst 
J 3 



from (^). Finally, I^^' rs should be reduced to the form (see e.g. fl44|)) 



(79) 



(80) 



r [d+],rs 
'4,* 



— I 

"PT 

\rsJt 



irs 
,0rs. 



jrs 
1 4 



E 

t=i 



irs 
Ms , 



jrst 
J 3 



(81) 
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The factor ^ in (|79| ) is the only occurrence of the inverse Gram determinant. (™) 
in (|80|) and fl81|) are in principle inverse Gram determinants as well, but since one does not 
expect the l.h.s. of ( |30D and ( |81~D to have kinematic singularities for (™J — > 0, the numerators 
also vanish in this case, which allows to keep the numerical evaluation under control. The 
situation is different in (f79| ) : even if the l.h.s. has no kinematic singularity as \£j — * 0, on 

the r.h.s. l| has been cancelled and therefore the above argument is not applicable. 
It is convenient to add up all contributions in fl78| ) by introducing the quantity 



W2 







With this definition (|78|) can be written in the form 

fOk\ 

\0r)t ^ 



n-JI [d+] 



E 

r=l 



6 

E 

s = l 
s^r 



{r} 



sO jrt 

ij M 



EW 



t=i 



st jrst 
ij 1 3 



(82) 



(83) 



As we see, some further cancellations occur. The representations derived here are particularly 
useful for numerical evaluation. 
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